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A closed-form solution for the potential in a current-free semiconductor 
surrounding a semi-infinite flat plate, carrying a small potential, is 
derived. We show that, jttst as in the case of the classical potential problem, 
there is a weak (square-root) singularity in the electric field at the edge of 
the plate. The solution also supplies an estimate of the effect of edge cur- 
vature, on edge field for a thick plate, found from equipotential contours. 
For a plate along the positive x-axis the normalized potential u is given by 



u = e 



erf c ft + y) + e - Jj _ erfc ft - g)1 



where erfc is the complementary error function; £, rj parabolic coordinates, 
such that ft + it}) 2 = x -f- iy; and x and y are in units of Debye length. 

I. INTRODUCTION 

In this paper we calculate the electrostatic potential in a current- 
free semiconductor surrounding a semi-infinite flat plate carrying a 
small potential. Our objective is to determine the nature of the sin- 
gularity in the field at the plate edge. Such information should be 
useful in the application of finite difference methods to similar, but 
more complicated, problems. We find that the field near the edge has 
the same weakly singular behavior as in the classical potential problem, 
being inversely proportional to the square root of the distance from 
the edge. 

It turns out, more or less fortuitously, that the present boundary 
value problem has a closed form solution, in terms of exponentials 
and error functions, so that it has some intrinsic mathematical interest. 
The solution was originally obtained by a very tortuous path. The 
problem was first attacked by the Wiener-Hopf technique. Then it 
was recognized that the inverse of the Fourier transform of the x- 
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derivative of the solution had a very simple form. After several changes 
of integration variable, integration of this derivative led to the closed 
form solution. 

Here we bypass this involved procedure and derive the solution 
directly, first factoring a particular solution and then transforming 
to parabolic coordinates. This after-the-fact derivation is admittedly 
artificial; in particular, up to now no other nontrivial boundary value 
problem has been solved by this "method." 

In the sections which follow we first derive the solution, then record 
various special forms, and finally consider its behavior in the neighbor- 
hood of the plate edge. The method of solution is discussed in Sections 
II and III; the results are discussed in Sections IV and V. 

II. THE BOUNDARY VALUE PROBLEM 

The potential u(x, y) due to a semi-infinite flat plate electrode in a 
current-free semiconductor or a quasi-neutral, stationary plasma sat- 
isfies the equation 

V 2 u = u xx -f w„„ = u, (1) 

the boundary condition 

u(x, 0) = 1, for x ^ 0, (2) 

for unit potential on the plate, the symmetry condition 

u v (x, 0) = 0, for x < 0, (3) 

and the limiting condition 

u(x, y) -> 0, (4) 

far from the plate. In the above the Debye length \ D has been taken 
as the unit of length and the plate potential — ip Q {< 0) as the unit of 
potential. The potential <p(X, Y) in dimensional form is then given by 

<p(X, Y) = -<p u(X/\ D , Y/\ D ). 

Equation (1) only holds for small electrode potential. Specifically, 
for an n-type semiconductor, we require that 

qcpo/kT « 1, 

for electronic charge q, Boltzmann constant k, and absolute tem- 
perature T. In this case the Debye length is given by 

for semiconductor dielectric constant c and donor number density N d . 
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III. THE SOLUTION 



Without preamble we now write down the solution to the above 
boundary value problem, justifying its form after the fact. Its structure 
is best displayed in a mixed notation. We find 

u = •"/({ + v) + e-'[l - /(| - ,)], (5) 

where £ and 77 are parabolic coordinates, such that 

x = £ - v 2 , y = 2^77, 
and the function / is given by 

2/(s) = erfcs = 1 - erf s = 2t~* I e~'' dt. 

By using the identities 

£ ± r, = r^cos- ± sin-j , 

we may also express u in terms of polar coordinates r = (x 2 + 2/ 2 )*, 
6 = tan -1 (y/x). 

Now let us attempt to unravel this rather involved expression for u. 
The occurrence of parabolic coordinates is not surprising. They are 
natural coordinates for flat plate problems with the family of parabolas 
77 = constant degenerating into the positive z-axis for 77 = and the 
orthogonal family £ = constant giving the negative z-axis for £ = 0. 
The way in which £ and 77 appear in equation (5) and the factors e ±v 
are a little more unusual. As we shall see, these two features are tied 
together. 

The exponential factors explicitly display the behavior of the po- 
tential far out along the plate, for fixed y and large x, where the potential 
should approach the one-dimensional potential e~ lwl . For our purposes 
they are also essential in simplifying the form of the differential equa- 
tion in parabolic coordinates. If we transformed equation (1) directly 
to parabolic coordinates, the resulting equation would be complicated 
by the presence of the "magnification factor" of the transformation, 
a function of £ and 77. On the other hand, with an exponential factor 
removed, the resulting differential equation does not contain the 
magnification factor. We have a special case of the following general 
result for equation (1): 

Suppose that u*(x, y) is any particular solution of equation (1). 
Then u(x, y) = u*(x, y)v{x, y) also satisfies equation (1) if v satisfies 
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the equation 

V 2 v + 2V(lnw*)-Vv = 0. 

Now, because the same factor occurs in both the Laplacian and in 
the product of gradients, this equation has the same form in all or- 
thogonal coordinate systems; that is, 

D lt + v„ + 2[(ln u*)^ + (In «*),(>,] = 0. 

In general, this identity is of no help, since In u* is a complicated 
function of £ and 77; but, for u* = e v , In u* = y = 2fr. Then we have 

»U + V„ + 4(7?^ + £y,) ■ 0. 
Notice the way that £ and 77 occur in this equation. If we set 

w(t, v) = /(* + l)> 

the equation reduces to the ordinary differential equation 
f (« + v) + 2(1 + *)/'(£ + 1?) = 0, 

an equation satisfied by /(£+ii) = constant and by /(£+*?) = erfc (£+17). 
A similar result obtains for the multiplier of e~ v , except that now we 
find a function of £ — 77, g(£ — 77), say. To satisfy the boundary condi- 
tion on the plate, we must have gr(£) = 1 - /(£). We are then left 
with two arbitrary integration constants in / which are determined 
from the conditions for fixed y and large positive and negative x. The 
result is equation (5). 

IV. THE POTENTIAL AND FIELD IN POLAR COORDINATES 

In polar coordinates equation (6) becomes 

e 
u = - 



(0) 



- — <1 — erf rM cos 5 + Bin 5, 

erf rM cos- — sin -J >• 



+ C — T— 11 + 



Various special values of u are of particular interest. Along the x-axis 
we have 

fl, for x ^ 0(6 = 0,r = x), 

u(x, 0) = «J 

erfc I x |*, for x ^ 0(0 = ir, r = \ x |), 
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while on the ?/-axis (0 = ir/2, r = y) 

m(0, y) = | erfc (2y)* + ^- 

Since for large y, erfc (2?/)* ~ (2iry)~*e~ 2v , this vanishes exponentially 
for large y, as it should. 

Near the edge of the plate (that is, for small r) 

u«l - 2(rA)*sin|- 

The a; and y components of electric field are given by 

u x = (irr)~ i e~ T sin - , (7) 

a 

u u = —(Trr)~*e~ r cos- 

i_ c+r8in9 /, f[i( *_l • e \]\ 

H ^ — l ? V C0S ^ sm W I 

- ^p^ {l + erf [r*(cos | - sin |)]}- (8) 

Note the simple form for w x , a separable solution of equation (1) in 
polar coordinates. The explicit form for u was in fact originally derived 
by first calculating u x and then integrating with respect to x, a rather 
laborious procedure. According to equations (7) and (8), the field near 
the edge of the plate is weakly singular; that is, | \7u | ~ (irr) for 
small r, just as in the case of the classical potential problem. Along 
the x-axis 



u x (x, 0) = 



u u (x, 0) = 



0, for x > 0, 

V \x \)~ h e~ [l \ for x < 0, 

— (irx)~*e~ x — erf (x) 7 , for x > 0, 

.0, for x < 0. 



Note that u u tends very rapidly along the plate to the one-dimensional 
value u v = — 1. The equipotentials near the plate then rapidly become 
straight and parallel to the plate. Any one of them can be chosen 
as an electrode in the form of a thick plate with the electric field large, 
but not infinite, at its edge, which now has a finite radius of curvature. 
Thus, at least for this family of plates, we can calculate the maximum 
field strength as a function of edge curvature. 
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V. THE EFFECT OF EDGE CURVATURE 

Suppose that x = x(y) is the equation of the equipotential u = u ; 
that is, 

u[x(y), y] = u„ . 
Then 

U x x' + M„ = 0, 

and 

u x x" + u xx x' 2 + 2u xu x' + u vu = 0. 

At the edge y = 

x = X (Q) = Xo , z'(0) = 0, 

where erfc | x |* = u . Thus the radius of curvature of the edge is 
given by 

R = 1/| z"(0) | = | u x (x , 0)/u uv (x o , 0) | 

= | u x (x , 0) |/| u(x , 0) - u xx (x , 0) |, 

using the differential equation (1). Now 

u(x, 0) = erfc | x |*, 

for x ^ 0, so that 

u x (x, 0) = e^'/Hx!)*, 

and 

w«(x, 0) = (1 +%\x\)u x (x,0). 

Thus 

12 - 2 |*. |/| 1 + 2 | *. |(1- 1/uOI, (9) 

where the normaUzed gradient 

W = I Vu(x , 0) \/u(x , 0) = e"' x,l /0r | a:„ |)* erfc | x |*. (10) 

If equations (9) and (10) are evaluated for various x , one obtains 
the edge field strength due to unit potential applied to a plate with 
the shape u = u , as a function of the edge radius of curvature, shown 
in Fig. 1. As 12 increases, one might expect that the edge field should 
approach the field on a circular conductor of radius R. This field, 
shown as a dashed curve in Fig. 1, is given by K 1 (R)/K (R), where 
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Fig. 1 — The edge gradient as a function of edge radius. 




Fig. 2 — Equipotentials around a flat plate in a semiconductor. 
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K and K x are the modified Bessel functions of the second kind. Except 
for very small edge radius, the fields very nearly coincide. 

Figure 1 also shows a sketch of a typical equipotential, that is, 
a thick plate of the family u = u„ . The shape shown is typical of 
moderately small edge radius R. For large R the plate thickness d 
tends to 2R, so that thick plates in the family tend to have semicircular 
noses. On the other hand, thin plates (u„ — > 1) have "wedge-shaped" 
noses, R tending to zero like (1 — u ) 2 , while d tends to zero like 1 — u„ . 

Figure 2 shows the exact shapes, for 0.1 ^ u„ ^ 1. Note the sharp- 
ness of the equipotential u = 0.9, compared with u = 0.1, even allow- 
ing for the difference in horizontal and vertical scales. For small edge 
radius then, unless the plate under consideration has a shape similar 
to the above, the present analysis is not likely to furnish useful estimates 
of the edge field. 

VI. ACKNOWLEDGMENT 

The author wishes to acknowledge many fruitful discussions with 
J. McKenna and E. Wasserstrom and, in particular, the computation 
of Fig. 2 by Miss Armida Marucchi. 



